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Abstract 


This report summarizes all the activity under the 
subject grant during the period 1977 December 1 - 1979 
May 31. The accomplishments during the period of the grant 
were to fit a first order analytical theory of Enceladus- 
Dlone to observational data, to develop a complete theory 
for Mlmas-Tethys , to develop new software for the calculation 
of the disturbing function In the case of Titan and Hyperion, 
and to develop the disturbing function Itself for Tltar. and 
Hyperion. 


APPENDIX I 


THEORY OP MIMAS AND TETHYS 


William )!• Jefferys and Lynne H« Rlea 
University of Texas at Austin 
Austin, Texas 78712, U.S«A< 


ABSTRACT 


An analytical theory of Mimas and Tethys has been 
developed, retaining terms which produce perturbations of 
the order of +10 km« The theory uses a novel set of 
variables, and was developed with the Hori-*Lie algorithm, 
using the algebraic manipulation language TRIGMAN. the 
perturbations have been implemented by means of FORTRAN 
subroutines produced by the computer^ 
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I. Accomplishments Under this Grant. 

a) Theory of Enceladus-Dlone 

Some 1,428 observations of Enceladus and Clone were fitted 
to the first order theory of Enceladus and Clone which was 
developed previously under NASA Grant NSG 44-012-282. The data 
consist of many right ascension and declination measurements 
on the two objects, as well as Information on orbital parameters 
reported by Kozal In Annals of the Tokyo Astronomical Observatory 
Series 2, Vol. 5, no. 2 (1957). The theory was extended to the 
second order of the small parameters, due to the fact that a 
number of slgnlf‘‘cant terms, both periodic and secular, exist 
In the second order which are Important at accuraclea which 
could be observed from space. 

b) Mlmas-Tethys 

The thecry of Mimas and Tethys was completed, and since the 
last semi-annual status report, a paper on this theory has been 
submitted to the Astronomical Journal. The theory of Mimas and 
Tethys represented a substantial advance over the Enceladus and 
Clone theory developed earlier. A new algorithm was developed 
for the calculation of the disturbing function for Mimas and 
Tethys, which was not only more accurate, but also consumed less 
computer time than the algorithm used for Enceladus and Clone. 

In addition, a new set of variables was developed, which had 
significant advantages over the variables used In the Enceladus 
and Clone Investigation. Indeed, we went back and recalculated 
the Enceladus-Clone theory using the new variables and found them 
to be advantageous In that case as well. Since the Mlmas-Tethys 
case Involves a large amplitude llbratlon. It was necessary to 
develop the theory to the 16th order In the small parameter 
of the llbratlon. This Is to my knowledge the first time that 
a 16th order development using the Hori-Lie algorithm has 
actually been applied to a practical problem. A copy of the 
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paper on Mimas and Tethys, which describes the theory in detail 
is attached as Appendix I of this report. 

c) Titan-Hyperion 

Since the last report, we have been concentrating on the 
theory of Titan and Hyperion. This theory has presented some 
unexpected new challenges, which have tested the TRICMAN 
algebratic manipulation system to its limits. Indeed, it 
turned out to be impossible with TRIGMAN alone to develop 
the entire disturbing function to the accuracy that would be 
required for a practical theory. The reason for this is that 
in the development of the disturbing function it is necessary 
to go to quite high order in the eccentricity of Hyperion 
in order to obtain a sufficiently accurate value of the motion 
of longitude of pericenter of Hyperion. Just to obtain accuracies 
comparable to the observed value, it is necessary to go to about 
the fourteenth or fifteenth power in the eccentricity. The 
development of terms in the disturbing function with such high 
powers of the eccentricity involves the subtraction of large 
quantities to produce small quantities, with consequent roundoff 
error. It turned out that about the tenth power of the eccen- 
tricity, the significance of the 60 bit coefficient word in 
TRIGMAN was completely lost. It therefore became necessary 
for a different strategy to be employed for the development of 
the disturbing function. 

The short period and Intermediate period terms in the 
disturbing function were developed as usual using the software 
which was used for Encealdus-Dione and Mimas-Tethys . The 
resonant and secular terms, on the other hand, were developed 
using entirely new software which has been developed during the 
past six months for this purpose. This involves calculating 
the particular terms that are required from a numerical Fourier 
analysis of special values of the disturbing function. In 
addition, certain derivatives of the disturbing function needed 
to be computed in the same way. These terms are then added to 
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the short and Intennedlate period disturbing function to 
produce the final disturbing function. In addition, since 
solar terms are of importance, they are added to the expression 
for the disturbing function. 

We now have in hand a definitive expression for 
the disturbing function, which is to be used for the theory of 
Titan and Hyperion. It is our intention to complete this theory 
in the near future, and to publish the results. We plan to do 
this on an unfunded basis. 

II . Publications 

One paper has appeared during the tenure of this grant. 

It is titled "Theories of Resonant Satellite Pairs in Saturn's 
System" by W. H. Jefferys and L. M. Ries, and was presented 
at I.A.U. Symposium y/81 in Tokyo, Japan May 23 - 26, 1978. 

A second paper has been submitted for publication, and is 
attached to this report. A paper on the application of the 
observations of Enceladus and Dione to the theory is in 
preparation, and one on the theory of Titan and Hyperion 
will be prepared after the calculations are complete. Copies 
of these papers will be made available to NASA as soon as they 
have been prepared. 

III. Staff 

During the lifetime of this grant, the staff consisted of 
W. H. Jefferys (Project Director) 3.15 man-months (1.25 man 
months paid from grant funds, 1.9 man months contributed time); 

L. Ries, 7.75 man months; N. Otto, 1.0 man months. During the 
final reporting period W. H. Jefferys contributed 0.6 man months, 
and L. Ries was paid for 2.0 man months. 
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1. INTR ODUCTIO N 

The near future pruDlscs greatly inproved precision In 
the observations of the satellites of the outer planets* 
Jupiter and Saturn flybys will produce greatly Improved 
data, and the j^O.002 arcsecund astrometric precision of 
Space Telescope, due to be launched In 1983, will also 
provide data of unprecedented accuracy* The 40*002 
arcsecond figure corresponds to an error of about 10 km at 
Saturn's distance* 

A program of regular observations of the satellites of 
the outer planets has been underway at McDonald Observatory 
since 1973 (Mulholland, Shclus and Abbott 1976, Mulholland 
and Shelus 1977, Benedict, Shclus and Mulholland 1978, 
Mulholland, Benedict and Shelus 1979)* At the same time, 
work on Improved orbital theories of Enceladus-Dlone, 
Mlnms-Tethys and Titan-Hypcrlon was commenced by the present 
autliorb (Jefferys and Rles 1973; Jefferys 1976; Jefferys and 
Rlcs 1979; hereinafter denoted Papers I, II and III, 

I 

respectively)* It has been the goal of these theories to 
provide terms to somewhat greater precision than that 
expected from the Improved observations which will soon 
become available* 

These theories are also of interest from a mathematical 
point of view* All possess a deeply resonant term In 
llbratlon, and tliere are also shallowly resonant terms* In 
this paper we describe the theory of Mimas and Tethys* 


”«go 3 


Since there arc aome important differences between this 
theory and the theory of Cnceladus and Dlone which was 
described In Papers I and II, It will be described in some 
detail > 

The theory of Mlmaa and Tethys, like the theory of 
Enceladus and Dlone, has been computed with the aid of the 
automated algebraic manipulation system TRIGMAN (Jefferys 
1970, 1972), using the Horl-Lle algorithm (Horl 1966), as 
formulated by Campbell and Jefferys (1970), This version of 
the Horl-Lle method Is particularly suited for automated 
computation. In order to carry the theory out to a 
sufficiently high order, however, we have found it necessary 
to Introduce a new set of variables, which has some 
Interesting properties. In addition, our theory has some 
other unique features In the way the deep resonance Is 
handled. 


II, VARIABLES 


In the theory of Enceladus and Dlone, we employed 
variables Aj, E^, 1^ (In the notation of Paper 1) that are 
linear functions of the usual elliptic elements (semlmajor 
axis, eccentricity and sine inclination). These had the 
advantage of simplicity, but led to a rapid increase In the 
size of series whenever Poisson brackets were calculated. 


This was due to the fact that the partial derivatives of the 
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variables Aj, Ej and Ij , with respect to the Dclsunsy 
variables, had to be expanded In powers of Aj, Ej and ^ 
themselves. Thli caused s rapid proliferation of terms In 
even the lowest orders, leading to s problem of 
"Intermediate swell" which made It difficult to obtain the 
complete solution. 


In the present work, we have used modified variables 
which have the property that the required partial 
derivatives consist of only a single term. The Delaunay 
variables are (Paper I) 




( 1 ) 


where J > 1,2 for the Inner and outer satellites, 
respectively. Dropping the subscripts for convenience In 


what follows, we have 
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where 

• 0 

Is the 

"nominal" value of the semlmajor axis 

and 

L 0 *^0 c> ^ 

• These equations have the 

solution 


a* - 

2i^ (L/L 0 -1) 


c* - 

(2(L - C)/I^ ) ^ 

(M 

1* - 

[2(G - H)/lo J ^ 



In order to write down the Hamiltonian we must express 
a, e and i in terms of the starred variables* This leads 
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Co: 


1, • Lo(l + / /2a o) 

2 2 

a - L /(w ) 

e - «j*<Lo/L^ (1 - 

2 k 

C - L(l - c ) ’ 

jL 2 L 

1 - l*(Lo/C) ’ (I - 1* 1^/4G) ^ 
H - C(1 - 1 ^ ^ 


(5) 


For the calculation, we employ variables A*, E^, and 1’*' 
which are related to a'’*^, e* and 1^ through the equations 


i* ■ a 0 + f 


e* - f F* 
e igfc 

1* - fjl* 


( 6 ) 


In a manner similar Co that used In Paper I. The 
Hamiltonian and all other quantities are expanded In powers 
of A*, E* and 1'*, and the truncation parameters 
and are used, as In Paper I, to control truncation of 

small terms* 


Our experience with 

these modified variables 

has 

been 

very gratifying* 

By 

cutting 

down on 

the 

amount of 

intermediate swell. 

the 

required 

computer 

time 

and 

field 


length have been considerably reduced, and the final results 
are quite manageable. 
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III. THE DISTURBING FU NCTION 

Our nat'od of calculating the dlaturblng function haa 
alao been Modified froa the Method ueed In Paper I. In that 
work, an Iterative achene ualng the Newton-Kapheon Method 
was used to develop the quantity 1 /a (the reciprocal of 
the distance between the two satellites). In the present 
work, It was decided that an aporoach using Lsplace 
coefficients %rauld be superior. Developing r| , rj and 
cos S as In Paper I, we write 

• 2 2 2 

‘A ■ rj + r2 - 2rjr2Cos S 

2 2 2 

Ag- aj + a2 - 2aja2C08 (A2 “Ai) (7) 

2 2 

A - Ao « 

and then expand 1/ A In powers of A , using Laplace 
coefficients to obtain the required powers of I/Aq. We 
have found that this method Is not only faster than our 
original method, but also more accurate. 


IV. NONRESONANT PERTURBATIONS 

There are some surprisingly large nonresonant 
perturbations In tne theory of Mimas and Tethys which appear 
not to have been computed before. For example, the 
argument ( 1 j + 2g j + 2h j - 21j - 2g 2 ” 2h ^ , which Is 
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responsijle for Clic deep rcsonanc* In titc cane of Kneel adtis 
and Dlone, haa a I'^afflclenl of over 8 ninutea of arc In tho 
theory of Hlmaa and Tethya, vltlch produces perturbations In 
the position of Mimaa of over 400 km. In addition, there 
are none saaller perturbations in twice this arRuuent, and 
perturbations in 1| which are due to the ohlatcness. To 
be sure, the principal perturbations in the Mlmas-Tethys 
theory, due to the deep resonance, are much larger, but 
these other terms are large enough to be observed, possibly 
even from the ground* 

In our work on Enceladus and Dionc, it was pointed out 
that TRICMAN, wliich is a Poisson scries processor, is unable 
to handle literal divisors, which necessitated a fairly 
complex scheme for cx;)anding the divisors in powers of small 
'i'.ntitics. This expansion was necessary, as was pointed 
out in Paper I, in order to calculate all of the 
perturbations correctly* In this paper we have adopted a 
similar approach, but have streamlined it and made it more 
systematic* The final results of our new approach are the 
same as if we had employed the method of Paper I, but the 
method is easier to carry out, and has the additional 
advantage that it postpones to higher orders the terms which 
contain high powers of tlie expansion parameters* 

The method we have employed in the present work is to 
adopt as the zero order Hamiltonian a function which, when 
differentiated with respect to the Delaunay variables, 
yields the numeric al values of th? corresponding mean 
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■otiona, rather than a aeries expancion* The original 
Haalltonlan for the problem is then partltluned (according 
to powers of the small quantities) among first, second, 
third, and higher order terms, and the xero-order 
Hamiltonian subtracted from the first order. Thus, If the 
original Hamiltonian of the problem Is 
F “ ^0*** ^2'*’ ***' * Hamiltonian %^lch 
reproduces the numerical values of the oiean motions, then 
the modified Hamiltonian that we use is given by 


Order 0: Fq - Fq 

Order 1: F'j - F^j+ Fj - Fq 

Order 2; F'2 - F2 (8) 

Order 3: F'3 - F3 


Slncu the modified Hamiltonian Is actually equal to the 
origlnui Hamiltonian when all the terms are added U|. , 
differing only In th-j way It Is partitioned among the 
various orders, the two systems roust have the same solution. 
However, for the original Hamiltonian, the equation for the 
elimination of the kth order periodic terms by the Hori-Lle 


method Cakes the form 



Pane 


wticre 


9W 

(Fo,V^> - E n^_V 




(9) 


(In the notation of Paper I), \A\ero the mean motlona nj 
are f unct Ions ; while for the modified Hamiltonian, the 
equation becomes 



3W' (10) 

where (F'.WM - E n i ^ 

Ok .1 — 


where the n j arc the nominal values of the mean motions. 
As a result, with the modified equationa, the divisors of 
the trigonometric terms arc always nu mbers . wlilch the 
Poisson series processor can handle, rather than 
expressions. 


The only drawback to this scheme is that in order to 
obtain the Solution to n given order (compared to the 
solution of the original Hamiltonian), it is necessary to 
carry the solution using the modified Hamiltonian one order 
higher. This turns out to be a small price to pay, as the 
series grow more slowly, and ’'cv terms are introduced Into 
the determining function only as tliey are needed. The 


algorithm for performing a Horl-I.lc transformation to 
arbitrarily high order is actually quite short when the 
method of Campbell and Jefferyr (1970) Is used - it takes 
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only about a half a page of TRIGMAN code* 

I 


t 


V. RESONANT PERTURBATIONS 

A'cer the elimination of the nonresonant part of the 
Hamiltonian, the method described In Paper II is used to 
handle the resonant perturbations. The basic approach is 
outlined In Paper III, and we develop It here in greater 
detail • 


The critical argument In the Mlmas-Tethys case Is 

21 + 2g + 3h - 41 - 4g ^ - 3lj , which Is replaced by an 

1 1 1 2 2 2 

angular argument, say 6 • Conjugate to It, a new action 
variable, 0 , Is Introduced by the method of Paper II, The 

Hamiltonian Is then expanded In powers of 0 to a 
sufficiently high order - In practice, only a few terms are 
actually needed, as 0 Is very small. The resulting 

Hamiltonian has the form 


2 

F - (AO - 2 B cos 0)/2 + OgO + ,,, 

2 2 

- (A 0 + B 0 ) /2 + Og 0 + , , , , 


( 11 ) 


where A and B are functions of the other variables and n„ , 

the mean motion of the critical argument. Is Identically 

zero. The condition n • 0 Is In fact a constraint which 

6 


must be satisfied In the final solution. It is one of two 
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constraints which must be Imposed to account for the two 
additional variables, 6 and e , which have been 
Introduced to handle the resonance. (The other constraint 
Is equivalent to the condition that the mean value of 8 
must equal scro.) 

In the second of Eqs. (11), wo have shown the coclne 
expanded In powers of 6 , and the lowest order term 
retained. In the theory of Enceladus and Dlone, where the 
amplitude of the llbratlon Is small, such an expansion Is 
quite adequate, and the hamlltonlan becomes that of a 
harmonic oscillator. In the theory of Mimas and Tethys, 
this Is no longer true. The amplitude of the llbratlon la 
very large, and many terms roust be taken In the expansion of 
the cosine in order to attain sufficient accuracy. 

The usual method of treating the Mlmas-Tcthys resonance 
Is, In fact, to use elliptic functions. However, TKIGMAN is 
not able to handle elliptic functions, and we have chosen a 
different approach, that of treating the problem as a 
perturbed harmonic oscillator (that Is, perturbed by the 
higher order terms In the expansion of the cosine). 

We replace 0 and 0 by a now canonically conjugate 
pair of variables p and q through the equations 

0- (2py)^ cos q. 

0 ■ (2p/ y) ^ sin q. 


( 12 ) 
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_ _ V 

where y “ (D/A)^ , and the bars Indicate Chat the functlona 
A and B have been evaluated at the noalnal values of their 
parameters. The Hamiltonian now takes the form 


h - — k 

F - p[A(B/A) + B(A/B)’ 1/2 

+ n (2py)^ cos q (13) 

+ plA(B/^^ - B(OT)^ 1 (cos 2q)/2 
3/2 

■f terms of order p and hlgner. 


Notice that the first term, at the nominal values of A 

— h 

and B, -is equal to (AB) , which is the nominal value of 

the mcai motion of the libration for small amplitudes. If 

it Is evaluated away from the nominal values, the mean 

motion will, of course be different. In fact, the 

coefficient of p In the first term of Equation (13) consists 
of the leading terras of the expansion of (AB)^ , which Is 

the actual mean motion of the small-aroplltudc llbratlng 
argument about Its nominal value (AB^^ . The second term 

of Equation (13) Is the one from which the principal part of 
the libration in longitude arises. As explained in Paper 
III, although the factor n^ vanishes because of the 
constraint condition discussed above, its partial 
derivatives do not vanish, so that this term, and any term 
factoted by a single power of , will give rise to 


nonvanishing perturbations In the angular variables. 
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The third tern of Equation (13) also la factored by a 
term that vanlahea, at least on the nonlnal orbit* This 
term can also give rise to perturbations In the angular 
variables, although they are much smaller than those arising 
from the second term. 

The Hamiltonian must now be put Into a form which Is 
convenient for the elimination of the angular argument q by 
means of the Horl-Lle method* As with the nonresonant 
perturbations, we take for the zero-order Hamiltonian a term 
which yields the numer i cal value of the mean motion of the 
llbratlon argument* The Hamiltonian Is then partitioned 
among first, second, and higher orders as before, using the 
power of the variable p to determine the order of a term* 
The Horl-Lle algorithm Is then applied, and the argument q 
eliminated from the Hamiltonian to the desired order* 

Because of the amplitude of the llbratlon, we have 
found It necessary to carry the solution out to the 

k 

sixteenth power of (p)^ * However, this should not be 

considered quite the same as a sixteenth order theory* What 
we are actually doing is nothing more than expanding the 
elliptic functions in trigonometric scries. The resulting 
determining function Includes terms through the argument 7q, 
and powers of p up Ihrougli the eighth* 

VT . PF.RTHKBATIOHS 1^ CO O RDI NATES 

In order to apply the theory. It Is necessary to 
compute the perturbations In tlie observed coordinates* With 
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the Hurl-L'..e method, this is very easy to do for any 
arbitrary function of the variables. The resonant and 
nonresonant determining functions are expressed in terms of 
transformed variables corresponding to the action-type 
variables A ^ and p, and the angular variables 1 , 

g , h and q , from which the short and long period terms 
have been eliminated. The transformed action-type 

variables, together with the zero points of the transformed 
angular variables, are the constants of Integration of the 
theory. They are close to the original variables E* and 
I* In terms of which the Hamiltonian was originally 

expressed. 

To obtain perturbations in any other quantity, such as 
the radius vector cr the longitude. It is only ncccessary to 
express the desired quantity In terms of the variables A*^, 
and , and the angular variables 1 , g and h . One 
then applies to this expression the Hori-Lle algorithm with 
the nnnresonant determining function to obtain the 

nonresonnnt perturbations. This function is then expanded 
in powers of 0 , wtiich in turn is expressed in terms of p 

and q by means of Equation (12). (The latter operation 

generally adds few, if any terms to the series, since 0 is 
so small.) 

The final step is to take into account the effect of 
the resonance by applying the Horl-hle algorithm to this 
scries, using tlie resonant (Ictcrmlning function. The result 
is a scries in the twice-averaged variables (which are 
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conacanta and linear functlona of the time)* Thla operation 
can be carried out for any dealred quantity which might be 
obaerved* We have done It for auch obaervablca aa the 
radlua vector, mean and true longitude, and latitude, aa 
well aa quantltlea auch aa the aemlmajor axla. Inclination 
and eccentricity* Perturbations In virtually any dealred 
quantity can be obtained eaally and efficiently with our 
programa* In addition, partial derlvatlvea In any quantity, 
with reapect to the conatanta of the theory, are caay to 
provide* 

It ahould be noticed that It la poaslble for “mixed" 
terma Involving both the llbratlon argument q aa well aa 
other angular arguments to appear In the solution* We find 
such terms In the present theory, and some of them are of 
significant size* For example. In the longitude of Mimas 
there Is a term with an amplitude of about 3 minutes of ->rc 
Involving both q and the shallow resonance 

The theory is available In the form of FORTRAN 
subroutines, automatically written by the computer, which 
are capable of calculating tlie perturbations In any of the 
functions mentioned above, such as longitude, latitude, 
radius vector, etc< Subroutines to calculate the partial 
derivatives of each of these functions with respect to the 
constants of integration arc also provided for the purposes 
of differential correction* It would not be difficult to 
provide such subroutines for any other appropriate function 


I 
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if that turna out to be necccasary* 
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